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We present the basic aspects of deep inelastic phenomena in the framework of 
the QCD parton model. After recalling briefly the standard kinematics, we discuss 
the physical interpretation of unpolarized and polarized structure functions in terms 
of parton distributions together with several important sum rules. We also make a 
rapid survey of the experimental situation together with phenomenological tests, in 
particular concerning various QCD predictions. In the summary, we try to identify 
Ch . some significant open questions, to clarify where we stand and to see what we expect 

to learn in the future. 
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I. INTRODUCTION 



\ The deep inelastic phenomena has been extensively studied, for about a quarter of a century, 

' both theoretically and experimentally. The main purposes of this vast physics programme were, 

\ first to elucidate the internal structure of the nucleon and later on, to test perturbative Quantum 

Chromodynamics (QCD) which generalizes the parton model. At the end of the sixties, the first 
' measurements on charged lepton deep inelastic scattering (DIS) at the Stanford Linear Accelerator 

,__( , (SLAC) have shown that the nucleon is made of hard point-like objects Q], which was the first 

■ evidence for the existence of sub-nuclear particles called partons |^ . Unlike the ep elastic scattering 

cross section which drops very rapidly for large momentum transfer, the observed cross sections 
for eN — > e' X were much larger than expected, when the off-shell photon probing the nucleon 
had a large . Actually these cross sections were independent and they obeyed the scaling 
Oh! behavior predicted by Bjorken ||^, whose physical picture was first given by Feynman in terms 

of the quark parton model (QPM). The explanation of scaling was based on the fact that a hard 
collision among partons takes place in such a short time, so that partons behave as if they were free 
objects. This significant observation led to the consideration of non-Abelian gauge field theories, 
which possess the crucial property of asymptotic freedom, and to propose QCD as the fundamental 
theory for strong interactions Q. It is now well established that QCD is an asymptotically free 
' theory and the strong interaction coupling UsiQ"^) becomes small when is large, i.e. at short 

distances. However scale invariance is broken because of quantum corrections and therefore one 
expects scaling violations which are calculable in perturbative QCD. 

The outline of this chapter is as follows. In section || we discuss the kinematics of DIS and define 
the relevant structure functions for both unpolarized and polarized cases. We provide the physical 



interpretation of the structure functions in terms of quark parton distributions in section III. 
We also give their main features from some recent data on unpolarized and polarized structure 
functions and we recall several important sum rules. In section ^ we go briefly one step beyond 
the QPM to consider scaling violations and to make phenomenological tests. Finally in section 
we give a short summary on significant open questions, and we try to clarify where we stand and 
what we expect to learn in the future. 
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II. KINEMATICS OF DIS AND STRUCTURE FUNCTIONS 



A. ep ^ e'X unpolarized case 

Let us start to recall the basic kinematics entering in the calculation of the cross section for the 
DIS process ep e'X to lowest order in Quantum Electrodynamics (QED), as shown in Fig. 
We define the Lorentz invariants in the usual way 

s^{p + kf, q^^{k-k'f^-Q\Q^>0), 

v^pq/M, =p'^ = + 2Mi^-Q^, (1) 

where M is the proton mass and we neglect the lepton mass. In the laboratory frame, if E and E' 
are the energies of the incoming and outgoing leptons and 9 is the scattering angle one has 

Q2 ^ 4:EE' sin^ 61/2 and v ^ E - E' , (2) 

which are the invariant mass and the energy of the exchange off-shell photon. One also defines two 
dimensionless scaling variables 

x^Q'^l2p.q = Q^l2Mv and y = p.q/p.k = v / E, (3) 

which can vary between zero and one. We recall that x is called the Bjorken variable and we note 
that X — 1 corresponds to elastic ep scattering, since from eq.(|l]), if — one has = 2Mv. 
Moreover for x = 0, i.e. = 0, we have real photoproduction. Let us recall that DIS corresponds 
to the ideal situation where x is fixed and and u are going to infinity, which is never achieved 
in practice. 




FIG. 1. The basic Feynman diagram for deep inelastic electron proton scattering in lowest order QED. 
The final hadron state which is not measured is denoted by X. 
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The unpolarized cross section described by the Feynman diagram shown in Fig. |^ has the 
following expression 



_ 4a^ E' 
dE'dfle ^ E 



^ III/ " J 



(4) 



where a = 1/137.036 is the fine-structure constant, L^/^} is the symmetric leptonic tensor 



(5) 



and where W^,y is the symmetric hadronic tensor defined as 



(6) 



X 



Because the electromagnetic current is conserved, namely 9^ = 0, we have q^W^u — quW^^i, = 
0, and by using parity conservation and T invariance, one shows that W^j/ can be expressed in 
terms of two real structure functions Wi and W2 



M2 



(7) 
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FIG. 2. World compilation of R as a function of for various small x values, curves are a fit of E143 
Collaboration [HI . 
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If one neglects the electron mass rUe, the cross section eq.(||) reads 
and for the clastic contribution when ~ 2Mv, i.e. a; = 1, we have with r = Q"^ /AM^ 



(8) 



tGIi5{v - Q^/2M) , 



tG 



M 



1 



Siu - Q^/2M) 



(9) 



where Ge and Gj\/ are the electric and magnetic form factors. Clearly this contribution vanishes 
in the large region since Ge,m{Q'^) ^ Q However for a point-like particle, from eq.(|^) 

we have 



2MT¥f * = x5{l - x) and vW^' = ^(1 - x) 



pt 



(10) 



where x — /2Adv, and we will come back later to this important result. From the 1968 SLAG 
experiment |jl|, it was observed that at a fixed value of W'^, the invariant square mass of X (see 
eq.(0)), when is large enough, both MWi and vW2 vanish, like in the elastic case. However in 
the scaling limit, that is when both v and are large, with x — Q'^ /2Mv fixed, this is no longer 
the case and according to Bjorken one expects 



MWi(v, Q2) 



Fiix) 



F2{x) 



(11) 



where Fi^2{x) are two scaling functions independent of Q^. Therefore the SLAG experiment has 
shown that the DIS cross section is much larger than expected and obeys scaling as predicted by 
Bjorken. These observations were the basis of the Feynman's idea for the parton model which can 
be formulated as follows: hadrons possess a granular structure and the "granules" behave as hard 
point-like almost free objects. 

Finally we recall that the structure functions Wi and W2 can be related to the absorption cross 
sections for transverse and longitudinal virtual photons ax and ctl 



W2{vM^) 
This allows to write eq.(8) as 



1 



:((TT + CTl) ■ 



(12) 



dE'dfle 

where s and F are kinematic factors 

1 



■tan^ 9/2 



27r2Q2 E 1-e 



and 



(13) 



(14) 



'^This behavior is firmly established for Gm{Q^) but recent Jefferson Lab data |^ have cast some doubts 
on its validity for Ge{Q'^) ■ 
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The positivity of (Tl, fT leads to the inequaUties 

0<Wi< (^l + ^^W2 . (15) 

One also defines a useful quantity, the ratio R of longitudinal to transverse cross sections for 
polarized virtual photons on an unpolarized target 

In the limit of large Q^, aL vanishes, so i? = and we have 2xMWi ~ vW2, that is the so-called 
Callan-Gross relation |^ 

2xF^(x) ^ F2[x) , (17) 

which is approximately valid for finite . However at finite Q^, quark transverse momentum 
effects generate non-zero values of R. These have been observed in the measurements made at 
SLAG showing that R is slowly varying with at fixed x. A recent world data compilation 
is depicted in Fig. g. 

Although the two scaling functions -^1,2(2^) were introduced above in the limit —> 00, for finite 
we have Fi^2{x,Q^) which arc dependent because of scaling violations. These important 
aspects will be discussed in section 0.1.4, but meanwhile in order to complete our presentation 
of the data, we give in Fig. || a compilation of Ff^i^iQ^) for fixed proton target experiments 
(NMC, BCDMS, E665) and collider experiments at HERA-DESY (ZEUS, HI). Within this broad 
kinematic domain in x and , the data exhibit indeed, a clear evidence for scaling violations. 
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FIG. 3. The structure function F'f{x,Q^), data are from refs. For the purpose of plotting, a 

constant c{x) — 0.6{ix — 0.4) is added to Ff 2{x,Q^), where ix is the number of the x bin ranging from 
ix = l {x = 0.32) to ix =21 {x = 0.000032) (Courtesy of R. Voss). 
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B. The neutrino induced reactions p(i')p /i^X 

For these reactions the proton is now probed by a weak charged current and the lowest order 
diagram is also shown in Fig. where 7* is replaced by a virtual W^. Actually, if are replaced 
by e^, one can also consider the line reversed reactions e^p v{p)X which will be studied at 
HERA DESY with polarized incident beams e^. The cross section is in this case 

(fa Gl ( Ml, V E' 



^ ^Li^M'\ (18) 



dE'dn {2tt)^ + J E 

with Gf — lO^^/A/^, My/ — 80.43GeV and where the leptonic tensor L^t'' is no longer symmetric 
Z/^'^-' — k^ki, + k^kjj — k.k g^^ it is^^pcrk^'k " , (19) 

with + for the left-handed neutrino and — for the right-handed antineutrino. 

For the hadronic tensor associated to the weak charged current, which is not conserved, one 
has now in addition to the symmetric part analogous to the previous case, an antisymmetric term 
due to the interference between vector and axial parts of the current. The hadronic tensor can be 
expressed in terms of three structure functions 

W^{l^,Q^) f p.q \f p.q 



M2 yP^-^i^)\P^-^i'^ 

WpVW^3''(^,Q') • (20) 



M 



The new structure function , which is not necessarily positive, is due to the parity violating 
coupling of the weak gauge bosons W^. The cross sections for v and 9 now read 



dE'dfl 27r2(l + Q2/M2,)2 

cos^ 9/2 + 2W'( SIT? e/2 T [ ] ^3 sin^ 



(21) 



and clearly by taking the difference between and P cross sections one can isolate W^. As in the 
previous case, one can interpret the W^{v, Q^) in terms of the virtual VF-proton cross sections and 
from positivity one obtains the inequalities 



< ^^^^J m\ < W( < (1 + u^lQ^)W^ . (22) 
Finally one also expects in the scaling limit 

vW^{v,Q^) F^{x), (23) 

where 3(2;) are three scaling functions independent of and similarly to eg. (|l7|) , one should 
have 
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2xFi'{x)^F,^{x)^x\F^{x)\ 



(24) 



For finite these three functions are also dependent and for completness we show in Fig. |4[ 
the results on xF^'^ from the CCFR collaboration @ at FNAL, as a function of Q^, for fixed x 
values. The existence of scaling violations is also shown by these data. 
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FIG. 4. The structure function 2;_F3'''^(a;, Q^), where A*' stands for nucleon, versus for different x 
values. Data are from ref. and for the purpose of plotting, a constant c{x) = 0.6(19 — ix) is added to 
xF^^ {x, Q^), where is the number of the x bin ranging from ix = 1 {x — 0.0075) to ix = 18 (x = 0.75). 



C. ep — > e' X polarized case 

Let us consider the case of DIS where both electron and proton are polarized. If Se denotes the 
polarization vector of the incident electron, it is a space-time vector orthogonal to the momentum 
k of the electron, so it is such that fc.Sg = and Sg — —1. Consequently the leptonic tensor 
has now, in addition to the symmetric part (see eq.(^)), an antisymmetric part 

= -2ie^,^0s:q^ . (25) 

Similarly, if Sp denotes the polarization vector of the proton the antisymmetric part of the hadronic 
tensor is linear in Sr, and reads 
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+ ^9" [{p.l)sl - Mp'] G,{., Q') , (26) 

where Gi^2 sue two spin-dependent new structure functions. These structure functions can be 
related to the absorption cross sections CTi/2 and (T3/2 for virtual photons with projection 1/2 and 
3/2 of the total spin along the incident photon direction, ct introduced above (see eq.(p^) is the 
total transverse cross section ctt = 5 ('^1/2 + '^3/2) and axL is the cross section corresponding to 
the interference between the transverse and longitudinal polarizations of the photon. We have 



lyMGi - 0^02 = — ^~2~~(^i/2 - o'3/2) and 



2 



02(AfGi + UG2) = „ aTL ■ (27) 



27ra2 

One can also define asymmetries for the virtual photon scattering, namely, 

0-1/2 - 0-3/2 vMGi - Q^G2 , 

Ai = — ■ '— = — and 

0-1/2 + 0-3/2 14^1 

= ^^^^ ^ y/Q^iMG^ + ,.G2) ^ ^28) 

0'l/2 + 0-3/2 Wl 

with the following positivity bounds [p^ 

|yli|<l and \A2\ < , (29) 



where R is defined in eq.(|16|). The first inequality is rather trivial unlike the second one, which 
can badly restrict the allowed values of A2. This bound has been recently improved as follows 



1^2! < V^(l + Ai)/2. (30) 

It involves Ai and in the kinematic regions where Ai is small or negative, one gets a much stronger 
bound than eq.(|29|) which corresponds to = 1. In the scaling limit one also expects 

iyM^Gi{iy,Q^) — > gi{x) and 

iy^MG2iiy,Q^) ,92(2:) , (31) 

where 31,2(2^) are two scaling functions independent of Q^. Of course in this limit, A2 is going to 
vanish since R also goes to zero. 

Finally these asymmetries Ai and A2 can be simply related to the measured spin asymmetries. 
In the case of longitudinal beam and target polarization along the beam direction (i.e. Se = 
— {\k\, Ek/\k\) and idem for Sp), one has 



where D and i] arc kinematic factors (the factor D is also called the depolarization factor) 

1 - E'e/E 



8 



In the region where v and are large, r\ is small and one has 

In the case of longitudinally polarized beam and transversely polarized target (i.e. Sp 
s.p = and s ^ = —1), one has 

where d and ^ are kinematic factors 



When V and are large 



^1 - dA, 



(34) 

(35) 

(36) 
(37) 
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FIG. 5. The world data on the asymmetries A2 for proton and deuteron targets taken from ref. [ p^ . 
The weak positivity bound ^TR is shown as the dotted curve. The solid curve was obtained by using 32 
given by the Wandzura-Wilczek sum rule, (see below eq.(|65|)) 

We show in FigJ^ a compilation of the world data on the asymmetries A2 for polarized proton 
and deuteron targets, which turn out to be rather small. Experimental results on the asymmetry 
A\ allow to extract the function gi, which is also dependent, and we show the most recent data 
at = 5GeV^ , for a polarized proton target in Fig. |6| and for polarized deuterium and neutron 
targets in Fig. 0. 



Although such a neutron target does not exist, the data are obtained from deuterium or from helium 3. 
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III. PHYSICAL INTERPRETATION OF STRUCTURE FUNCTIONS, PARTON 
DISTRIBUTIONS AND SUM RULES 



A. The physical picture : the quark parton model 

The experimental facts observed at SLAC in 1968 ||l[ gave a strong support to the existence of 
point-hke objects in the proton structure and were the inspiration for trying to elaborate a simple 
physical picture for interpreting the structure functions of DIS. When E oo (i.e. v — s- oo), 
because of time dilatation, the life time of the proton virtual states is also very large. When 
oo, the interaction time 1/\/Q^ becomes very small. So in the scaling limit, v and very 
large, the virtual photon 7"^ (see Fig. 0) sees a "frozen" state of quasi-real, quasi- free point like 
objects called partons. 




FIG. 8. The parton model handbag diagram showing the interaction of 7* with the parton of momentum 

k. 

According to the parton model, in this kinematical region, the proton behaves as a gas of free 
partons. In particular, it is very unlikely to recover a proton in the final state and this is why the 
elastic form factors Ge,m vanish when 00, as recalled above. So in the calculation of the 

DIS cross section {e.g. see eci.(||)), the hadronic tensor W^^ is evaluated according to the parton 
model graph, shown in Fig. y, namely 

^/^^ = E E / '^^''fsiP' k)S[ik + q)% (38) 

i s •' 

where fl{p, k) is the probability for the virtual photon to struck a parton i of spin s and momentum 
k, inside the proton of momentum p. It seems very natural nowdays to assume that the charged 
partons which interact directly with 7* are quarks^. This step was not so obvious to make in 
1968, although we had already the SU{3) classification of hadrons, as built up from quarks, and 
the Cabibbo theory of weak interactions simply formulated in terms of leptons and quarks. It is 



^Gluons are also fundamental constituents of the proton, but since they don't carry electric charge, they 
don't couple directly to 7*. 
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important to recall that nowdays six quarks have been identified. They have very different masses: 
three of them are the light quarks up, down and strange (u, d, s) and three are the heavy quarks 
charm, bottom and top (c, 5, t). The quark flavor corresponds to one of these six labels u, d, s, c, b, t 
and the u, c, t quarks have a charge e = 2/3, while the charge of the d, s, b quarks is e = —1/3. The 
summation over i in eq.(^8|)) and subsequent equations stands for a flavor summation i = u,d, s, ... 
In the nucleon u and d quarks dominate and their mass can be neglected in high energy DIS. So 
if these partons are indeed spin 1/2 massless quarks, we are now in the framework of the quark 
parton model (QPM), ^/^^^^{q, k) is the tensor which describes the interaction of 7* with the quark 
i of momentum k and we have, for its symmetric part 



e2 



w 



(q,k)^^Tr{^j^{^ + 4)j,} , (39) 



which is similar to the leptonic tensor {e.g. see eq.(^). If is the fraction of the proton momentum 
carried by the quark i, Xi = k.q/Mv, so SYJt+qY] — j^S{x~Xi) and the hadronic tensor becomes 

i s 

[2kf,k^ + kf,q^ + k^q^- g^i^k.q], (40) 

where is the charge of the quark i. From this expression, one can extract VFi,2 whose features 
correspond to the scattering of point- like objects (see eq.(|lO|)) and by comparing with eq.(0), one 
finds after some partial integration and using eq.(|l l|) 

^i(^) = ^E^'*(^) ' (41) 

i 

where 

l^(^) = E / d'f'Tflip, k) = q,+ {x) + (42) 

Clearly qi{x), which is called the unpolarized quark distribution function, contains the sum over two 
quark spin directions and the average over the quark intrinsic transverse momentum fc^. Usually 
in the QPM one interprets qi^(^_^{x)dx as the probability to find a quark i, with helicity parallel 
(antiparallel) to the nucleon spin, and with momentum fraction of the proton between x and x-\-dx. 
From eq.(pO[) one also finds VW2 = 2xMWi, that is the Callan-Gross relation (see eq.(p7|)) which 
leads to i? = 0, as mentioned before. This is what one expects for quarks of spin 1/2 and let us 
recall the argument. Consider a collision quark-virtual photon in the Breit frame. By conservation 
of the total spin along the momentum direction, a scalar quark (s = 0) cannot absorb a photon 
of helicity A = ±1, so ctt = and R — 00. However for s — 1/2, the quark helicity flips in the 
electromagnetic interaction and therefore only the photon with A = ±1 can be absorbed, so ctl = 
and i? = 0. 

Finally let us consider polarized DIS and the interpretation of the polarized structure functions 
51,2(2;). Similarly in the QPM, the antisymmetric part of the hadronic tensor (see eq.(p6|)) can be 
expressed as 

W^M^ -EE '^'''fsiP' 'pHt^'(l^ k, s)S[ik + q)% (43) 
i s 

where Sp denotes the proton polarization vector and 

wli^' = ^Tr{{l + + 4h^} ■ (44) 
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So 

W^i^^ = i^t.i'paqp d'^k [flip, k, sp) - fi {p, k, Sp)] ■ 

i 

k„5[{k + qf] . (45) 
By comparing with eq.(26) and using eq.(|3l|) one finds 

9i{x) = i^e,?Ag,(x) and 32(2:) = , (46) 

i 

where with the help of the definitions in eq.(^) 

/\q,{x) ^ qi+{x) - q,^{x) , (47) 

which is caUed the polarized quark distribution function. Therefore in the scahng hmit the longi- 
tudinal polarization of the proton is described by gi{x) and we find that 32(2;) vanishes. 

Actually if one goes beyond the most naive QPM (see section |l^) , one must also consider the 
existence of gluons which produce quark-antiquark pairs. Since the antiquark contribution can 
be computed similarly, all the above quark distributions must be replaced by adding quarks and 
antiquarks, so the correct expressions are indeed 

Pi{^) (2;) + (2;)] and 

i 

91 (2^) = ^ E (^) + (^)] • (48) 

i 

One sees from these expressions that by measuring only Fi{x) (or i^2(2)), it will not be possible to 
disentangle neither unpolarized quark distributions of different flavors, nor quarks from antiquarks. 
For polarized quark (antiquark) distributions the situation is the same if one measures only gi{x). 
So we give now a short discussion to clarify this point. 

Let us go back to y{p)p — > p^ X which are interpreted in the QPM in terms of four different 
scattering processes, i.e. vq, vq, vq and Dq. For vq scattering {resp. uq) both u and q are left- 
handed [resp. both v and q are right-handed), so the cross section da/dy is isotropic (s wave) 

whereas for vq since q is right-handed, in this case {p wave) da/dy ~ (i±f^) - (1 - ^J)^ {tdem 
for i?q since q is left-handed). Therefore one can write for a given flavor 



dxdy 
dxdy 



[qix) + (1 - yfq{x)] and 

[q{x) + il-yfq{x)] . (49) 



Consequently the y dependence allows to separate q{x) from q{x) and by measuring both v and 
P cross sections, one can isolate q{x) — q{x) and therefore we have (see eqs. (|l|) and (H)) 

xF^ix) = ^(«*(^) - ■ (50) 

i 

A similar analysis can be done to separate Aq{x) from Aq{x) by using charged current reactions 
in e^p polarized DIS, which might be done at HERA DESY in the future. 
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B. Main features of quark distributions from DIS data 



Once the parton structure of the proton has been detected, to estabhsh it firmly, we need to 
make a detailed study of its basic constituents. A very high energy proton is not a simple three 
quarks objects like a proton at rest, but much more complex and it is actually made of valence 
quarks, gluons and quark-antiquark pairs which carry a small fraction of its momentum (see Fig. 
^). We define the valence quark distributions as 

qval{x) = q{x) - q[x) = q(x) - Qseaix) , (51) 

where qsea{x) denotes the sea quark distribution and we assume qseai^) = qi^)- From neutrino- 
nucleon DIS via the charged current, one can extract valence quarks from the measurement of 
F^{x) (see eq.(|50|)) and the most accurate data comes from CCFR Q shown in Fig. |. By 
looking at a fixed value, as shown in Fig. |l^, we see from this pure valence 
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FIG. 9. A microscopic view of the proton seen by a 7* having a large and a large u. 

contribution, that xu^ai{x) and xdyaiix) are maximum around x ~ 0.2 or so and they vanish for 
a; '--^ 1 and x 0. Since one can also isolate the antiquark distributions, the results from CCFR 
p4| show that xq{x) (or xqseaix)) goes to zero very rapidly for x ^ 1 and has a sharp rise for cc ~ 
(see Fig. [Tl| ). The small x kinematic region which corresponds to fixed and — > 00, is called 
the Regge region. A simple Regge analysis of the absorption cross sections 7*p leads to 

qsea{x) — > l/x and qvai{x) — > l/%/5 . (52) 

Clearly the singular behavior shown in Fig.|ll| is stronger than this prediction and we will come 
back to this important point. From ep —> e'X and ed — + e'X DIS on hydrogen and deuterium 
targets, one can extract the electromagnetic structure functions -Ff (a;) and F2{x) for proton and 
neutron, in terms of valence quark and antiquark distributions. 



14 



1 



0.8 



X 0.6 



0.4 



0.2 



10 



10" 



FIG. 10. The structure function xF^'^ {x) versus a; at = 3GeT/^ taken from ref. 
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FIG. 11. The antiquark distribution at = 3GeV^ (full circles) and — 5GeV^ (full triangles) versus 
X taken from ref. |l4|. 
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By using eq.([l7|) and eq.(^8|) we have 



^ii^) = n [^"i^vaiix) + 2xu{x)] + - [xdyai{x) + 2xd{x)] 

+ -[xs{x)+xs{x)] , (53) 
9 

if we restrict ourselves to three flavors u,d and s. F2{x) is obtained from eq.(p3[) by exchanging u 
and d. The NMC experiment at CERN M provides a rather accurate determination of F2 (x) / F2 (x) 
and — F2{x) as shown in Fig.[l4 These data confirm and complete some of the features 

observed in neutrino data. The difference F2{x) — F2{x) has a maximum around a: = 0.3 or 
so, where Uyai{x) and d^aiix) dominate and it vanishes for a; ~ 1. As we will see later, its 
behavior in the small x region has revealed the flavor symmetry breaking of the sea quarks, i.e. 




FIG. 12. Ratio F^/F^ as a function of x for differents values, data are from NMC Q and E665 
[|o|. Difference F| - Fa" as a function of x for = 4GeV^ data are from NMC g. 

Since sea quarks dominate for small x, we expect the ratio to go to one in the 

low x region as indicated by the NMC data. The E665 data, which have larger errors, suggest a 
departure from this trend in the very small x region. The behavior of F2{x)/F2{x) for large x 
indicates that Uyai{x) dominates over d^aiix). The NuSea-experiment has measured the ratio 
of muon pair yields from Drell-Yan production of proton-proton or proton-deuteron interactions. 
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This ratio allows the determination of d/u ratio (see Fig. 13), which has a strong x dependence 
and confirms the antiquark flavor asymmetry, mentioned above and the fact that d > u. Using a 
parametrization of J + w, one can also extract d — u (see Fig.^4|), which is mainly positive (except 
may be at large x) and exhibits a rapid rise for low x, in agreement with what was observed in 
neutrino DIS (see Fig. [Tl] ). 




FIG. 13. d/u a.tQ = 7.35GeV/c measured by FNAL E866 experiment 



FNAL Ease 

= 7.35GeV/c 



FIG. 14. d-u &t Q = 7.35GeV/c extracted from FNAL E866 experiment [|l|. 
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The simple parametrization d{x) — u(x) = (1 — x)"^ , which was proposed ten years ago |^^, is 
consistent with these new data (see FigJIj) . The antiquark flavor asymmetry is an important topic 
which has generated new developments and for a complete review see ref. [ p3[ . 

We have already said that a direct flavor decomposition of the polarized quark distributions is 
not possible from the data on gi{x) and gf{x) (see eg. (^8|) ) , but this can be done using polarized 
semi-inclusive DIS, where one detects a hadron in the final state. We show in Fig. |l^ some recent 
results from HERMES for the determination of xAuv{x), which turns out to be positive, 
xAdv{x) which is negative and xAu{x) which is compatible with zero. The accuracy of these 
results is limited, but these trends agree with some global fits of the g^'"' {x, Q^) data, as we will 
see later. 

We now turn to a short review of different sum rules, the structure functions we have introduced 
so far must satisfy. We will also discuss their validity or breakdown versus experiment, as well as 
some theoretical considerations. 




0.03 0.1 0.5 



X 

FIG. 15. The polarized quark distributions at — 2.5GeV^ separately for the valence quarks xAu^{x), 
xAd^{x) and the sea quarks xAuv{x) as a function of x (Taken from ref. jsij). The solid lines indicate the 
positivity limit and the dashed lines are a model parametrization. 

C. Sum rules 

There exists a number of sum rules for unpolarized and polarized structure functions, some of 
which are rigorous results and other which rely on more or less well justified assumptions. Let us 
first consider the unpolarized case and for the electromagnetic structure functions we start with 
the Gottfried sum rule |p5| (G SR). If one assumes an SU{2) symmetric sea, i.e. u(x) — d{x), one 
can easily show, using eq.(p3h, that 
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Ia= f'^[Fiix)-F-ix)] = \, (54) 

Jo 'l^ -3 



because 



Nu ^ dx Uyai (x) = 2 and Nd = dx dvai (a;) = 1 . (55) 
Jo Jo 

In fact the NMC experiment [|j has observed a large defect of the GSR, since their measurement 
(see Fig. gives at = 4GeV^, Iq — 0.235 ± 0.026. This flavor symmetry breaking, more 
precisely d > u, is a consequence of the Pauli exclusion principle which favors dd pairs with respect 
to uu pairs, since the proton contains two u quarks and only one d quark. 

Next for the charged current structure functions, we have two rigourous results based on eq.(p5|), 
namely the Adler sum rule ^6|] (ASR) 



1 dx 
2x 



[Ff (x) - {x)]=N^-Nd^l (56) 



(an experimental verification | |37| with large errors show an agreement) and the Gross- Llewellyn 
Smith sum rule (GLSSR) 

Igls^ ^[xF;'P{x)+xF^''{x)\^N^ + Nd^?,. (57) 
Jo 

The ASR is exact and receives no QCD corrections, but its experimental verification is at a very 
low level of accuracy ||3^. The GLSSR gets a negative QCD correction and the CCFR data (see 
Fig.|) gives Igls = 2.55 ± 0.06 ± 0.1 at ^ 3GeV^ in fair agreement with the theoretical 
prediction. If we turn to polarized structure functions, there is first a fundamental result called 
the Bjorken sum rule (BSR). It was derived about thirty years ago in the framework of quark 
current algebra and it relates the first moment of the difference between gi{x) for the proton and 
gi{x) for the neutron and the neutron /3-decay axial coupling 

' dx[gP{x)~9'^{x)]^^gA/9v , (58) 



where gA/gv = 1.2573 ± 0.0028 is very accurately known. The BSR gets also a negative QCD 
correction and we will come back later to the test of this firm prediction of QCD. One can also 
derive sum rules for g^ and separately. These are the Ellis-Jaffe sum rules (EJSR) which 
read 

1 

T\= dxgl{x) ^ —{QF - D + QAs) and 
Jo 18 

r?^ = ^ dxgl (x) = ^ (6F ~ + 6 As) , (59) 

where F = 0.459 ± 0.008 and D = 0.798 ± 0.008 are the /?-decay axial coupling constants of the 
baryon octet and As = j}As{x)dx is the total polarization of the proton carried by the strange 
quarks. One recovers eq.(pq) by taking the difference because F + D — gA/gv- In their original 
work, Ellis and Jaffe made the critical assumption that As = 0, which allows to make definite 
predictions for and F". 



19 



^60 
a 

§40 



20 




- / \ 




a) 


















m 


















2 4 


6 , ,8 






[GeV^] 




6 2 28 

[GeVI 

FIG. 16. a) Igdh for proton, as a function of Q'^ for various upper limits of integration < 4.2GeV'^ 
(triangles), < 45GeV^ (squares), the total integral (circles) from HERMES col. M. The curve is a 
theoretical prediction b) IgdhQ^ /{iGi^^ce) as a function of shows the weak Q -dependence of Fi. 

Using the spin-dependent photoabsorption cross sections o'i/2(3/2){^)j Gerasimov-Drell-Hearn 
PI (GDH) have derived the following sum rule valid for real photons 



GDH 



[0-1/2 (i^) - 0-3/2 (i^)] 



1^' 



(60) 



where v is the photon energy in the target rest frame, i^thr is the pion production threshold and 
K the anomalous nucleon magnetic moment. This important sum rule which predicts — 204/x6 for 
the proton, has never been tested directly due to the need for a circularly polarized beam with a 
longitudinally polarized target and a wide range of photon energies to be covered. However the 
GDH integral can be generalized to the case of absorption of polarized transverse virtual photons 
with 



nOO 1 

IgDh{Q') = / ['Ji/2{'^,Q')-<J3/2{'^,Q')] — . 



One can show that, within a good approximation, one has in the scaling limit 



/GDff(Q') 



(61) 



(62) 



where Fi = jj^ dxgi{x), so the GDH sum rule is connected to polarized DIS. The Q^-dependence 
of the generalized GDH sum rule has been measured by the HERMES Collaboration B2I in the 
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range 1.2 <Q^< 12GeV^ (see Fig.|l6|). The Q^-behavior of Igdh{Q^) agrees with the theoretical 
prediction of ref. and this suggests that there are no important effects from either resonances 
or non-leading- twist. 

Concerning the structure function .92(2^), which is related to transverse polarization (see eq.(^7|), 
but has no simple interpretation in the parton model, it is possible to derive a superconvergence 
relation by considering the asymptotic behavior of a particular virtual Compton helicity amplitude. 
This leads to the Burkhardt-Cottingham sum rule P4| (BCSR) 



dxg2{x) — and 



dxg^ix) = , 



(63) 



for proton and neutron and from this result, it has been naively argued that 172(2;) vanishes identi- 
cally. Actually from what we discussed in section || about the asymmetry A2 in the scaling limit, 
one can alternatively expect gi (x) -\- g2{x) = 0, a simple relation between gi and g2. However g2{x) 
is more complicated than that Bs] and only part of it (twist-2 contribution) is entirely related to 
gi {x) by means of the Wandzura- Wilczek sum rule (WWSR) which reads for J > 1 



dxx 



.1-1 



J 



g^{x)+g^^{x) 







Clearly for J = 1 one recovers the BCSR eq. 



(|63| ) and for J = 2 one has 

9i{y)— - ■ 

V 



(64) 



(65) 



Very preliminary data obtained recently by E155X and shown in Fig.^ for proton, are consis- 
tent with eq.(p5|), as shown by the dashed curves, but this has to wait for further confirmation. This 
approximation gives also the same kind of agreement with the A2 data for proton and deuteron, 
as noticed above in Fig.|^. 
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FIG. 17. The structure function versus for different x values. Preliminary data from SLAG 
E155X H. 
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IV. BEYOND THE QPM, SCALING VIOLATIONS AND PHENOMENOLOGICAL 

TESTS 

A. Scaling violations 

So far we have considered the proton structure in the framework of the QPM, a picture vaHd 
only in the scahng hmit. For finite Q^, there are several experimental facts which confirm the 
existence of scaling violations and we have already shown in Figs. ^ and B the dependence of 
F|'(a;, Q^) and xFi^^ {x, Q^), for different values of x. For small x, ^^(a;, Q^) increases with larger 
and for large x, it decreases with larger Q^. Actually this can be understood in the framework 
of QCD, a color gauge theory, which is the underlying field theory for describing quark interactions 
and in particular, deep inelastic phenomena. 

From the knowledge of the structure functions , F2, F!^^ and F^", one basic observation, which 
was made long time ago, is the fact that in DIS, the total fraction of the proton momentum carried 
by quarks and antiquarks, i.e. j /o^ (^) is about 1/2. So half of the proton 

momentum is carried by electrically neutral objects, which have been identified with the gluons, 
already mentioned above. Gluons are the carriers of the color force and therefore in the QPM 
diagram shown in Fig.|| (see also Fig.^, one can no longer neglect the possibility that the quark 
may radiate a gluon, before or after being struck by the virtual photon 7*. These contributions 
are of order aas(Q^), where as(Q^) is the strong interaction coupling and moreover a gluon as a 
parton in the proton can contribute to DIS via the diagram 7*3 — s- qq to the same order. From 
the QCD Lagrangian and the corresponding Feynman rules one can compute these contributions, 
so partons in QCD are no longer described by distributions which are functions of x alone, 
but of X as well as Q. These scaling violations will affect the predictions of the QPM and, as an 
example, the ratio R (see eq.(|6|)) is not equal to zero, so the Callan-Gross relation (see eq.(pT|)) 
is not exact any more, and so on. The dependence of the parton (quark, antiquark and gluon) 
distributions has been investigated by several authors |49[ and can be summarized in the following 
way. Given the parton distributions at some reference point — Q^, the evolution is given 
by a set of coupled equations between singlet part^ q^ = q + q and gluon G 

d ( q^ix,Q^)\ ^ asiQ') 

f'dzf P^^{z) NfP,G{z) 
J, z \Pg,{z) Pgg{z) 

where Nf is the number of flavors and the Pij{z), the so-called splitting functions, are known 
[^,|9|. This set of equations, usually called the DGLAP equations, is valid in leading order (LO) 
QCD, but it can be generalized to next-to-leading order (NLO), i.e. by including corrections of 
order {as{Q^))^, because we have 

P,^ix)=pg\x) + ^P^^^{x) + .... (67) 

A more formal approach to derive the behavior of the structure functions is based on the use of 
the operator product expansion (OPE). The starting point is to consider time ordered products 
of two currents, which are related to the hadronic tensors introduced above, and to study their 
singularities on the light cone |5l[ | . This is also a very powerful method to derive directly sum rules 
in QCD, together with their radiative corrections at the level of one or several loops. 



*There is a similar equation for the non-singlet q^^ — q — q with no gluon contribution. 
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Concerning the polarized distributions there exists a set of equations similar to eq.(|6^) which 
provides their evolution [Q, and the corresponding polarized splitting functions can be found 
in refs. [^,^. We note that in addition to the quark (antiquark) polarizations Aqi(a;,Q^), one 
also introduces the gluon polarization AG(a;, Q^) defined as in eq . (p7|) . This distribution plays also 
a crucial role in connection with the anomaly of the axial current |34| , as we briefly recall now [^5| . 
For a given flavor i, the first moment Aq.; is defined as the nucleon matrix element of the quark 
axial current 

< p, s|gi7^75gi|p, s >= 2spAgi . (68) 

Due to the fact that this axial current is not conserved, it was shown that eq.(|6|) gets a gluon 
contribution which generates a shift of Aqi 

Aq, ^ Aq, - ^^f^AGiQ') . (69) 
Ztt 

It is worth recalling that from the QCD evolution equations Aqi is independent, whereas 
AG{Q^) increases logarithmically with Q^, but such an increase exactly compensates the decrease 
of as{Q^) with Q^. However the size of this gluonic correction, which can be partly absorbed in 
the definition of the quark distributions, is purely a matter of scheme convention [ p6| . Finally we 
note that so far, we have no direct experimental information on the sign and magnitude of AG. 



B. Phenomenological tests 

There is an enormous literature devoted to phenomenological work on DIS but we shall discuss 
only very briefly, a few aspects of this broad area of high energy physics. Clearly one particular goal 
is to parametrize all the parton distributions as function of x, at a given low value, = Qq, 
and to test whether the evolution agrees with the QCD predictions. In several analysis the 
quark and gluon parton densities were described by a polynomial standard form 

xp{x, Ql) = T^Ax^il - xf{l + 7a; + /5%/J) , (70) 

where t], a, 6, 7, p are free parameters. These twenty or so parameters were determined by making 
a global NLO fit of the data on unpolarized structure functions and the results obtained are rather 
satisfactory. 

In a different approach, one is using, for valence quarks and antiquarks, a set of Fermi-Dirac 
distributions | ]58[ |, with the following expression 

Ax"' 

xp{x,Ql) = , (71) 

exp[(x — Vyx))ix\ + 1 

where V{x) plays the role of the "thermodynamical potential" and x is the "temperature". The 
gluon density function G{x, Qq) is given in this approach by a Bose-Einstein distribution with no 
free parameter. 

In this case we have a small number of parameters, say ten at most, which were determined at 
= 4GeV^ by fitting the CCFR data displayed in Fig. |l^ for valence quarks, in Fig. [ill for 
antiquarks and the NMC data on the ratio and the difference of F2 and in Fig. 12, 
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From these distributions at = 4GeV^, after a standard evolution using the DGLAP 
equations at LO only H, one obtains the unpolarized parton distributions for different flavors, 
shown in Fig. 18 at Q — 20GeV^. We note that, as expected u quarks dominate over d quarks, 
d(x) > u{x) and the s and c quarks are very small. Clearly the gluon distribution G{x, Q^), which 
was divided by a factor 10 and represented by the curve g/10 in the Figure, dominates largely 
in the small x region. The predictions obtained in a larger range are in excellent agreement 
with the CCFR data for xF^^ {xjQ"^) depicted in Fig. and also with the very recent data at 
HERA-DESY from the Zeus and HI collaborations on FP{x,Q^) displayed in Fig. | The 

rapid rise in the very small x region, down to x = 10^^^ or so, where sea quarks dominate, is 
well described by the QCD predictions and this rise is the direct extrapolation of what has been 
observed by CCFR for x = 10^^. We gave only these two examples for illustration, but these tests 
have been widely extended as discussed in ref. [p8[. The large amount of available data on the 



unpolarized structure functions was also used to perform a global analysis by several authors |57 
and they all find for the parton densities a general pattern similar to the one displayed in Fig. 18 
We now turn to some tests involving the polarized structure functions. As we have seen before, 
there is a substantial experimental information on the x dependence of with proton, deu- 

terium and ^He polarized targets, for in a rather limited range, i.e. 

2 < Q2 < lOGeV^ Given 

the present accuracy of the data and the small range investigated so far, the evolution of 
gi'"{x,Q'^) predicted by QCD cannot be tested at the level it was done for the unpolarized case. 
Moreover this test depends on AG{x, Q^), for which only theoretical speculations can be made. In 
ref. [5^ ], since the initial quark (antiquark) distributions are qi±{x,Q^) for each flavor i, one gets 
immediately Aqi{x, Q^). Therefore one can calculate consistently the F2's and the gis structure 
functions, with no additional parameters and the agreement with the data is rather good [Q. We 
show in Fig. |l^ a typical pattern of polarized parton distributions. A.u(x) and Ad{x) are rather 



well determined and they are consistent with the results shown in Fig. [15[ We also see that A-O(x) 
is small and we find that Aq(x), for q — d, s,c are compatible with zero. The AG{x) shown in the 
Figure is one possible solution but it has large uncertainties. There are several other QCD analysis 
of the gi's p9| , where one relates the polarized parton distributions to the unpolarized ones with 
some additional parameters. These analysis conclude that AG is weakly constrained by present 
data and they all agree that a small AG is prefered. 

Finally let's say a few words on testing the first moment of gi, i.e. the EJSR and the BSR. 
Since data are obtained for each cc at a different Q^, in order to have gi(x, Q^) for a constant 



one has to rely on the crucial hypothesis that the asymmetry Ai (see eq.(28)) is independent, 
so the evolution of gi is driven by that of F2 . 

Moreover, the sum rules involve integrals from x = to x = 1 and since there is no data at 
both ends, one has to make some reliable extrapolations in these two regions, low x being the most 
critical one. The QPM results given in eqs.js^) and ( p9| ) are modified by QCD corrections, which 
have been calculated up to order {as{Q^)). For example, in the gauge invariant factorization 

scheme, perturbative QCD corrections to r^*-"-* lead to the following results |Q : 

r?^"^ = C^s (±^g\ + l^g\) + ^ Csg% (72) 



where 



g^^ AT, = Au + Ad + As, 

g\ = F + D, (73) 

g\^3F~ D, 



new work updating ref. |58| with new data and a NLO evolution, is in preparation. 
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the coefficients non-singlet Cns and singlet Cs are given by 



Cs 



as 




as\ 


TT 


12 V 


n / 




- 1.10 








TT 




V TT 



20.2153 



for three quarks flavors. The BSR becomes now 



1 5a 

Qgv 



TT 



12\ TT 



20.21531 — 

TT 



(74) 



(75) 



Here we have neglected the theoretical uncertainties in the very low region due to the higher- 
twist effects which go like 

The world data for the Fi (Q^)'s show a defect in the EJSR for proton and deuterium both at 
Q2 = 3GeV^ from E143 and = lOGeV^ from SMC |l|,||. From this one deduces also 

a defect for neutron, which was confirmed by the E154 result. If one combines the world data one 
finds for the BSR 



r^-r5' = 0.171 ± 0.005 ±0.01 at Q^ = 5GeY^, 



(76) 



to be compared to the QCD prediction 0.182 ± 0.005. One can also turn the argument around and 
given the experimental result, use the theoretical expression for the BSR to extract as{Q^) []6l| . 
We also note that the defect in F^ and F" are usually interpreted as being mainly due to a large 
and negative As (see eq.(59)), i.e. As — — 0.13±0.02, but of course this does not necessarily holds 
if one considers SU{S) symmetry breaking effects The total contribution of the quark spins 
to the proton spin, is defined as AS = Am + Ad + As. For As = 0, the Ellis- Jaffe prediction 
gives AS ~ 0.58, but according to all experiments, one gets a smaller value, for instance, E155 
experiment |27| has obtained AS = 0.23 ± 0.04 ± 0.6 at — 5GeV^, which is not so easy to 
reconcile with our first intuition from the naive quark model, namely, AS — 1. 



V. SUMMARY 



Our knowledge of the unpolarized and polarized structure functions has improved considerably 
over the last twenty years or so and consequently also, our basic understanding of the proton 
structure. The enormous quantity of high statistics experimental data which is now available, 
allows us to make rather precise tests of the QCD parton model and this theoretical picture is 
nowdays established on very serious basis. Deep inelastic phenomena is only one of the several 
areas of high energy particle physics which are the testing grounds of perturbative QCD and 
obviously, it has to be supplemented at least by electron-positron annihilation and hard hadron- 
hadron collisions [Q. 

For the unpolarized parton distributions, the flavor symmetry breaking of the sea quarks {i.e. 
dsea > Usea) was Only discovcred recently with the violation of the GSR and it should be in- 
vestigated more precisely for example in Drell-Yan lepton pair production or in production 
p8[ . Our knowledge about the gluon distribution, i.e. its small x behavior and dependence, 
is improving but a real experimental effort remains to be done to put it on equal footing with the 
quark distributions. It requires, for example, better measurements of -Fl* (x, Q^) in the very small 
X region and this makes HERA-DESY really exciting for the coming years. 

Concerning the polarized parton distributions significant progress have been made and we have 
now the following results. The violation of the EJSR is confirmed for the proton and for the 
neutron. The BSR is verified up to 10% and this important test must be strengthened at a 
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higher level of accuracy. The kinematic domain presently accessible is still too limited to answer 
some important questions, in particular about the small x behavior of gi, the evolution of the 
quark distributions and about the, so far, elusive AG{x, Q^). The perspective of having polarized 
protons at HERA-DESY, which is seriously envisaged ||6^ , is opening a new experimental window 
for answering these questions. Finally let us mention the very exciting spin physics programme 
which has been approved at RHIC-BNL and will start operating by 2001 as a polarized pp collider 
p5| . Hopefully it will also shade some light on parton helicity distributions as well as quark 
transversity distributions. A facility called eRHIC, which combines a few GeV polarized electron 
beam onto one of the high energy polarized proton beam, is also under serious studies at BNL. 



W. K. H. Panofsky 1968, rapporteur talk. Proceedings of the 14*^ International Conference on High 

Energy Physics, Vienna, p. 23 (Editors J. Prentki and J. Steinberger, CERN Geneva). 

J. D. Bjorken 1969 Phys. Rev. 179, 1547. 

R. P. Feynman 1969 Phys. Rev. Lett. 23, 1415. 

H. D. Politzer 1974 Phys. Reports C14, 129. 

A Collaboration, M. K. Jones et al. 2000 Phys. Rev. Lett. 84, 1398 and references therein. 
C. G. Callan and D. G. Gross 1969 Phys. Rev. Lett. 22, 156. 
L. W. Whitlow et al. 1990 Phys. Lett.B 250, 193. 
SLAG E143, K. Abe et al. 1999 Phys. Lett. B452, 194. 

New Muon GoUaboration, M. Arneodo et al. 1994 Phys. Rev. D50, Rl and references therein; 1995 
Phys. Lett. B364, 107. 

ZEUS Collaboration, M. Derrick et al. 1995 Z. Phys. C65, 379; 1996 ibidem C72, 399; J. Breitweg et 
al. 1997 Phys. Lett. B407, 432. 

HI Collaboration, T. Ahmed et al. 1995 Nucl. Phys. B439, 471; S. Aid et al. 1996 ibidem B470, 3; 
C. Adloff et al. 1997 ibidem B497 3. 

BCDMS Collaboration, A. C. Benvenuti et al. 1989 Phys. Lett. B223, 485. 
FNAL E665 Collaboration, M. R. Adams et al. 1996 Phys. Rev. D54, 3006. 

CCFR Collaboration, P. Z. Quintas et al. 1993 Phys. Rev. Lett. 71, 1307; W. C. Leung et al. 1993 

Phys. Lett. B317, 655; J. H. Kim et al. 1998 Phys. Rev. Lett. 81, 3595. 

M. G. Doncel and E. de Rafael 1971 Nuovo Cimento 4A, 363. 

J. Soffer and O. V. Teryaev 2000 Phys. Lett. B490, 106. 

SLAG E155 Collaboration, P.L. Anthony et al. 1999 Phys. Lett. B458, 529. 

EMC Collaboration, J. Ashman et al. 1989 Nucl. Phys. B328, 1. 

Spin Muon Collaboration, D. Adams et al. 1994 Phys. Lett. B329, 339; Erratum, 1994 Phys. Lett. 

B339, 332; 1997 ibidem B396, 338; B. Adeva et al. 1997 Phys. Lett. B412, 414. 

SLAG E80 Collaboration, M. J. Alguard et al. 1976 Phys. Lett. 37, 1261; 1978 ibidem 41, 70. 

SLAG E130 Collaboration, G. Baum et al. 1983 Phys. Rev. Lett. 51, 1135. 

SLAG E143, K. Abe et al. 1995 Phys. Rev. Lett. 74, 346. 

SLAG E143, K. Abe et al. 1995 Phys. Rev. Lett. 75, 25; 1998 Phys. Rev. D58, 112003. 
SLAG E142, P. L. Anthony et al. 1993 Phys. Rev. Lett. 71, 959; 1996 Phys. Rev. D54, 6620. 
Spin Muon Collaboration, D. Adams et al. 1995 Phys. Lett. B357, 248; 1997 Phys. Rev. D56, 5330; 
B. Adeva et al. 1998 Phys. Rev. D58, 112001. 

SLAG E154 Collaboration, K. Abe et al. 1997 Phys. Lett. B405, 180. 
SLAG E155 Collaboration, P. L. Antony et al. 2000 SLAG-PUB- 7994 ( |hep-ph/000724^ ). 



HERMES Collaboration, K. Ackerstaff et al. 1997 Phys. Lett. B404, 383; A. Airapotian et al. 1998 
Phys. Lett. B442, 484. 

P. Amaudruz et al. 1995 Nucl. Phys. B371, 3. 
D. Adams et al. 1995 Phys. Rev. Lett. 75, 1466. 



27 



[31] FNAL Nusea Collaboration, E. A. Hawker et al. 1998 Phys. Rev. Lett. 80, 3715; J. C. Peng et al. 1998 

Phys. Rev. D58, 092004. 
[32] G. Proparata, P.G. Ratcliffc and J. Soffer 1991 Phys. Rev. Lett. 66, 687. 
[33] S. Kumano 1998 Phys. Rep. 303, 183. 

[34] HERMES Collaboration, K. AckerstafT et al. 1999 Phys. Lett. B464, 123. 
[35] K. Gottfried 1967 Phys. Rev. Lett. 18, 1154. 
[36] S. L. Adler 1966 Phys. Rev. 143, 1144. 

[37] WA25 Collaboration, D. Allasia et al. 1985 Z. Phys C28 321. 
[38] D. Gross and C. Llewellyn Smith 1969 Nucl. Phys. B14, 337. 
[39] J. D. Bjorken 1966 Phys. Rev. 148, 1467 and 1970 ihid Dl, 1376. 
[40] J. Ellis and R. Jaffe 1974 Phys. Rev. D9, 1444; 1974 DIO 1669 (E). 

[41] S. B. Gerasimov 1966 Sov. J. Nucl. Phys. 2, 430; S. D. Drell and A. C. Hearn 1966 Phys. Rev. Lett. 
16, 908. 

[42] HERMES Collaboration, A. Airapetian et al. 2000 Phys. Lett. B494, 1. 

[43] J. Soffer and O. V. Teryaev 1995 Phys. Rev. D51, 25. 

[44] H. Burkhardt and W. N. Cottingham 1970 Ann. Phys. 56, 453. 

[45] X. Ji 1995 Contribution to the Workshop on Deep Inelastic Scattering and QCD, Paris April 24-28, p. 

435 (Editors ,J. F. Laportc and Y. Sirois). 
[46] S. Wandzura and F. Wilczek 1977 Phys. Lett. B 72, 195. 

[47] SLAC E155X Collaboration, contribution to the workshop on the Nucleon Structure in High x-Bjorken 

Region, P. Bostcd et al. 2000, Philadelphia (USA), March 2000. 
[48] G. AltarclU 1982 Phys. Rep. C81, 1. 

[49] V. N. Gribov and L. N. Lipatov 1972 Sov. J. Nucl. Phys. 15, 438,675; G. Altarelli and G. Parisi 1977 

Nucl. Phys. B126, 298; Yu. L. Dokshitzer 1977 Sov. J. JETP 46, 641. 
[50] K. G. Wilson 1969 Phys. Rev. 179, 1499. 
[51] R. A. Brandt and G. Proparata 1971 Nucl. Phys. B27, 541. 
[52] W. Vogelsang 1996 Nuc. Phys. B475, 47; 1996 Phys. Rev. D54, 2023. 
[53] R. Mcrtig and W. L. van Neerven 1996 Z. Phys C70 637 and references therein. 
[54] G. AltarelU and G. G. Ross 1988 Phys. Lett. B212, 391; R. D. Carlitz, J. C. Collins and A. H. 

Mueller 1988 Phys. Lett. 214, 229; A. V. Efremov and O. V. Teryaev 1988 JINR Report E2-88-287 

unpublished. 

[55] For a more complete discussion on this and on many other aspects of polarized DIS see M. Anselmino, 
A. V. Efremov and E. Leader, 1995 Phys. Rep. 261, 1; B. Lampe and E. Reya, 2000 Phys. Rep. 332, 
1. 

[56] For an extensive review on the issue of factorization schemes see H.Y. Cheng 1996, Int. J. Mod. Phys. 
A 11, 5109. 

[57] A. D. Martin, R. G. Roberts and W. J. Stirling 1994 Phys. Rev. D50, 6734; A. D. Martin, R. G. 
Roberts, W. J. StirUng and R. S. Thorne 1998 Eur. Phys. J. C4 463 and references therein; CTEQ 
Collaboration, J. Botts et al. 1993 Phys. Lett. B304, 159; H. L. Lai et al. 2000 Eur. Phys. J. C12, 
375; M. Gliick, E. Reya, M. Stratmann and W. Vogelsang, 1996 Phys. Rev. D53, 4775; T. Gehrmann 
and W. .J. Stirling 1996 Phys. Rev. D53, 6100. 

[58] C. Bourrely and J. Soffer 1995 Nucl. Phys. B445, 341. 

[59] T. Gehrmann and W. J. StirUng 1995 Z. Phys. C65, 461; R. D. Ball, S. Forte and G. Ridolfi 1995 
Nucl. Phys. B444, 287; F. Buccella, O. Pisanti, P. Santorelh and J. Soffer 1996 Nuovo Cim. A109 159; 
C. Bourrely, F. Buccella, O. Pisanti, P. SantoreUi and J. Soffer 1998 Prog. Theo. Phys. 99, 1017; D. de 
Florian, O. Sampayo and R. Sassot, 1998 Phys. Rev. D 57, 5803; E. Leader, A. V. Sidorov and D.B. 
Stamenov 1998 Phys. Rev. D 58, 114028; Y. Goto et al. 2000 Phys. Rev. D 62, 034017; M. Botje 2000 
Eur. Phys. J C 14 285. 

[60] S. A. Larin and J. A. M. Vermaseren 1991 Phys. Lett. B259, 345; S. A. Larin 1994 Phys. Lett. B334, 
192 and references therein. 

[61] J. Ellis and M. Karliner 1995 Phys. Lett. B341, 397; Workshop on the prospects of spin physics at 
HERA, Desy-Zeuthen, August 28-31, 1995 (Editors J. Bluemlein and W. D. Nowak) DESY-95-2000, 
p. 3; G. Altarelli, R. D. Ball, S. Forte and G. Ridolfi 1997 Nucl. Phys. B496, 337. 

[62] J. Lichtenstadt and H. J. Lipkin 1995 Phys. Lett. B353, 119; B. Ehrnsperger and A. Schafer 1995 



28 



Phys. Lett. B348, 619. 
[63] G. Sterman et al. 1995 Rev. Mod. Phys. 67, 157. 

[64] Proceedings of the Workshop on the prospects of spin physics at HERA, Desy-Zeuthen, August 28- 
31, 1995 (Editors J. Blucmloin and W. D. Nowak) DESY-95-2000. Proceedtngs of the Workshop on 
"Polarized Protons at High Energies - Accelerator Challenges and Physics Opportunities", DESY 
Hamburg, May 17-20, 1999 ( Editors A. De Roeck, D. Barber, G. Raedel). 

[65] N. Saito, invited talk at the 14th Int. Spin Physics Symposium, SPIN2000, Oct. 16-21, Osaka, Japan, 
to be published in the Proceedings. 



29 



